Parity breaking and scaling behavior in the spin-boson model 
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We study the breaking of parity in the spin-boson model and demonstrate unique scaling behavior 
of the magnetization and entanglement around the critical points for the parity breaking after 
suppressing the infrared divergence existing inherently in the spectral functions for Ohmic and sub- 
Ohmic dissipations. Our treatment is basically analytical and of generality for all types of the bath. 
We argue that the conventionally employed spectral function is not fully reasonable and the previous 
justification of quantum phase transition for localization needs to be more seriously reexamined. 
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The spin-boson model (SBM) has been key to phe- 
nomenological descriptions of open quantum systems, in 
which the environment acts as a bosonic bath responsible 
for dissipation of the system, i.e., the spin [J, Besides 
the coherence of the spin, the correlation between the 
spin and the bath degrees of freedom has also attracted 
much attention. 

The SBM hamiltonian is given by 0] 



H 



k 



^Afc(4-f-afc)CT^, (1) 



where CTz and ax arc usual Pauli operators, e and A are, 
respectively, the local field (also called c-number bias 0) 
and the tunneling regarding the two levels of the spin. 
a\ and ak are creation and annihilation operators of the 
bath modes with frequencies w^, and Afe is the coupling 
between the spin and the bath modes, which is governed 
by the spectral function J(aj) = X]fc '^^'^('-^ ~ ^k) for 
< a; < Wc with the cutoff energy Wc- In the infrared 
limit, i.e., u — >0, the power laws regarding J{uj) are of 
particular importance. Considering the low-energy de- 
tails of the spectrum, we have = 2T:au]].~''ijj'' with 
Q < u < LJc and the dissipation strength a. The expo- 
nent s is responsible for different bath with super-Ohmic 
bath s >1, Ohmic bath s =1 and sub-Ohmic bath s <1. 

There are several approaches solving the SBM, such 
as the non- interacting blip approximation 0, numeri- 
cal renormalization group (NRG) j34Tl|, quantum Monte 
Carlo (QMC) and so on [ij, Il4|. The main concern 
in the SBM is for the localization of the spin, and quan- 
tum phase transition (QPT) between the delocalization 
and localization in the Ohmic dissipation has been well 
investigated so far But the second-order QPT with 
sub-Ohmic dissipation, which is currently under intensive 
investigation, is not yet fully understood. 



In contrast to the intensively studied QPT, we inves- 
tigate the breaking of parity in the SBM. We show that 
variation of the parameters in Eq. (1) leads to differ- 
ent symmetries of the SBM hamiltonian and the parities 
to be broken are responsible, respectively, for localiza- 
tion and delocalization. The key step in our treatment is 
the suppression of the infrared divergence existing in the 
spectral functions for Ohmic and sub-Ohmic dissipations, 
which enables us to demonstrate unique scaling behavior 
of the magnetization and entanglement in the vicinity of 
critical points for the parity breaking. More importantly, 
our treatment is basically analytical and suitable for all 
types of the bath, by which we can fully understand the 
physics behind the infrared divergence and the scaling 
behavior. 

We start from suppression of the infrared divergence 
in the spectral functions. With reference to the standard 
form of the spectral function J(w), we introduce a dis- 
tribution function p{u}) = 1 — e^^("/"<^) with P a very 
large number, which is a smooth variation in the function 
of ijj with p(0) = 0. So the spectral function is modified 



J'(w) ^ TT^ A|(5(a; - Wfc)p(a;) = 27ratj^-''wV(w), (2) 



which fits J{lo) very well, as shown in Fig. 1 in the case of 
P — 10®. However, using J'{uj), the values of integration 
near the zero frequency could be effectively suppressed 
due to the exponential factor in p(w). 

To check how well the modified spectral function 
works, we compare calculations in the following with 
J{uj) and J'juj). We solve the SBM using displaced co- 
herent states [o, Ell as the eigenfunction of Eq. (1), 
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where c, 
later and {n} 



{„} and d{„} are coefficients to be determined 



ni, 



, n N are for different bosonic 
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FIG. 1: (color online) Comparison between J'{ijj) and J{uj) 
with P = 10^ J'{uj) fits J{u}) very well with the dilTerence 
hard to be distinguished from the curves. 



modes. is the product of displaced co- 

herent states of different modes, i.e., = 

nf=i l"fc>A,(B,)> where 



with the displacement variables qk = Xk/ijJk and k = 
1, 2, • • • , N. Using Schrodingcr equation, we have, in the 
case of A <c 1, 




|^^a;fe(TOfe - ql) + 



(3) 



\^u:u{mk - ql) - | 



A 



{m}j 

(4) 

where other terms, except D{,n,m}, in i^im.n} have been 
neglected due to the reasons in Supplementary Material 
m- ^{m,m} is given by [13, 111 



k=lj=0 



[{■nik-jWi^ 



It is straightforward to yield following solutions 
from Eqs. (3) and (4), that is, the eigenener- 

gics i?f„} - T.k^k{mk - ql) ± ^JT^^~^/2, 
and the coefficients c^m} ' 



and df„} = l/y^l + (/^f„,})2 with Aif„} 



e ± 



A2D2 



{rn,T7i} 



It is obvious from the ex- 



pression of eigenenergies that the ground-state energy 
E'^Q^ is smaller than E^^^. 



The above analytical treatments for Eq. (1) can be 
considered as complete and reliable solutions for the char- 
acteristic of the SBM [l^. For our purpose, we may fo- 
cus on the ground-state characteristic of the model to 
see counter-intuitive phenomena in the SBM. In such a 
case, the infrared divergence is reflected in the variable 
Z?{o,o}i which is written as i?{o,o} = e"^^*"^* [2^. We 
first check ql in the case of the bath modes of the con- 
tinuous spectrum 0, 0] ■ From the conventional spectral 
fmiction J{uj), we have 



— |<5(w — u)k)du} 



(5) 



with 



if.<l 

P={Hzr,) if^ = i 

if^>i> 

where a and Wc are defined above in the spectral function. 
uji is a small quantity regarding the frequency difference 
from u! = 0. In the case of the infrared limit, i.e., uji — >■ 0, 
we have J2kll °° ^ ^ which is actually caused 
by the uncertainty in the spectral function for uji — >0. 
However, if using the modified spectral function J'{uj) 
and repeating Eq. (5), we have ql = 2aj3' with 



1 

s-l 



s-l 
2 ' 



2P(s-i)/2 

ln(P) + r(o, P) 



P) 



if s = 1, 



where 7 is the Euler-Mascheroni constant, r(-) is the 
gamma function and r(-, •) is the upper incomplete 
gamma function. Since the gamma functions are finite 
even for a very large value of P, /?' is definitely con- 
vergent. Similar results could be obtained for the bath 



modes of the discretized spectrum [21 



We have noted the results in previous publications that 
there arc QPTs in Ohmic and sub-Ohmic dissipation 
cases, but not in super-Ohmic one, which exactly cor- 
responds to the infrared divergence demonstrated above: 
Divergence for Ohmic and sub-Ohmic bath, but conver- 
gence for super-Ohmic bath. As a result, it is reason- 
able to presume that the QPT presented previously are 
probably induced totally or partially by the infrared di- 
vergence in the calculations using J{lu). In fact, there 
have been some discussions about the shortcomings in 
NRG methods, which cause the qualitatively incorrect 
results when studying quantum-critical phenomena, and 
spoil the determination of critical exponents and behav- 
iors Q . Additionally, there were hints that the NRG dis- 
plays truncation errors or other errors in the long-range 
ordered phase H, [13, [H, [H . 

To fully understand the characteristic of the SBM, we 
may first consider two special cases of Eq. (1), where 
we denote the case of e 7^ with A = (e = with 
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IS an even 



A ^ 0) by F(e 0,A = 0) {H{e = 0,A 7^ 0)). For the 
two special hamiltonians, we introduce, respectively, two 
parity operators 11/ = az and 11// = o-^je"^^'' "fc"''. For 
H{e ^ 0,A = 0), we have [7J(e 7^ 0,A = 0),n/] = 0, 
with the ground state of their common eigcnfunction to 

be li^lo) = (j{0})^ satisfying H/ = - \^-^), i.e., 

an odd parity state of 11/. Since {iPYqIctz iV'Fo) ~ 
IV'Jo) always a localized state. Similarly, we have 
[H{e = 0, A 7^ 0),n//] = 0. The ground state of their 

common eigcnfunction = ^ (^|{0}j"* 

parity state of Ft// with Ujj jV'j^ q) = jV'j^ q). We have 

{''PTi ol '^^ \''l^7i 0) ~ meaning lipjj g) to be always a de- 
localized state. It is evident that the odd (even) parity 
breaks in the variation from A = with e 7^ (e = with 
A 7^ 0) to both A 7^ and e 7^ because the hamiltonian 
H in Eq. (1) never commutes with any of the parity oper- 
ators above. This also means that the ground state of H 
would never stay forever in delocalization or localization, 
but possibly moving between delocalization and local- 
ization in variation of certain characteristic parameters, 
such as a. We show below the behavior of magnetization 
in the vicinity of critical points of the parity breaking. 

The magnetization of the SBM is of importance to 
symbolize the transitions between delocalization and lo- 
calization of the model. Using the modified spectral func- 
tion and the displaced coherent states \{n}) a{B) , we have 



(6) 



where the average is made by the ground-state of the 
model and k ~ e/A. Since is convergent for any type 
of the bath, {<Tz) should be of finite value. Performing 
the second derivative of (cr^) with respect to a, we obtain 
a reflection point ac = — ln(2K^)/8/3', by which Eq. (6) 
is rewritten as 

(7) 



under the scaling transformation a' ~ a/uc- For a fixed 
value of K, the magnetization in Eq. (7) is only relevant 
to a', rather than other characteristic parameters. So etc 
can be regarded as a scale of the dissipation strength. In 
addition, if wc set a' = 1, the magnetization turns to be 
a constant —1/^/3, which implies a fixed crossing point 
for different types of the bath in the magnetization with 
variation of a' (See Fig. 2(a,b)). 

It is more interesting to demonstrate the scaling be- 
havior of the magnetization with a displaced dissipa- 
tion strength a" = {a - ac)l3'/V^- Since (crz(a")) = 



-1/v 1 4- 2e~24\/3a" ^ which is independent of both n 
and s under the scaling transformation, the magnetiza- 
tion with respect to a" , as presented in Fig. 2(c,d), re- 
mains unchanged for different types of the bath and dif- 
ferent tunneling and localization parameters. The scal- 
ing transformation was usually used to find QPT around 
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FIG. 2: (color online) The scaling of the magnetization, with 
(a): as a function of a' under sub-Ohmic dissipation for dif- 
ferent ft; (b): as a function of a' for different k and types 
of the bath; (c) and (d); as a function of a", which remain 
unchanged for the characteristic parameters k and s in the 
model. 



the critical points, where the scale invariance appears 
in the neighborhood of the critical points. In contrast, 
our results present the scale invariance in the whole re- 
gion of a", which can be understood as the critical be- 
havior resulting from the parity breaking regarding 11/ 
and 11//. It could be more clarified if we check the lin- 
ear variation near the region of a" = with the slope 
d{az) /da"\a"^o = —8 (See Fig. 3), which is a continu- 
ous change between the delocalization and the localiza- 
tion without any cusp-like behavior. 

It was indicated in previous studies for the Ohmic 
damping at e = that quantum Kostcrlitz-Thouless tran- 
sition separates the delocalizedphase at small a from the 
localized phase at large a @, [ll]. In contrast, the situa- 
tion of e = in our case only corresponds to delocaliza- 
tion and there is no possibility for any QPT. However, for 
Eq. (1), there is possibility of translation (with no cusp- 
like behavior) between the localization and delocalization 
in our results, where the delocalization and localization 
correspond, respectively, to small a {a < ac) and large 
a (a > ac). Nevertheless, our results is only relevant to 
the critical behavior of the parity breaking and hold for 
not only the sub-Ohmic damping but also other types of 
the bath. 

The scale ac has some unique features: It is universal 
for different types of the bath, which means a continuous 
and smooth curve with respect to s (See Fig. 4(a)). In 
contrast, if we employ the conventional spectral function, 
/3 in the expression of ac would yield ac ^ in the case 
of s <1, but finite values for s > 1, which causes drastic 
changes in the variation of ac with respect to s and corre- 
sponds to the appearance of QPT around the point s = 1 
(See Fig. 4(b)). This is another evidence that the QPT 
for the localization in the SBM is related to the infrared 
divergence. On the other hand, e = is a singularity in 
ac, as shown in Fig. 4(c), which is, as mentioned above. 
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FIG. 3: (color online) The magnetization in variance with a" 
in the nearby region of a" = under different characteristic 
parameters. 



Eq. (8) is replaced by \fTIo" - ln(2K2)/8 under the scal- 
ing transformation. Since the magnetization reaches for 
a < Q!c, i.e., the delocalization, and drops to -1 if a > ac, 
i.e., the localization (Refer to Fig. 2(c,d)), we could know 
from Fig. 5 that the delocalization and localization in the 
SBM correspond, respectively, to the increasing entangle- 
ment and the decreasing entanglement. In this sense, the 
scale ac is also the reflection point for the entanglement 
increasing and decreasing. As a result, we may easily con- 
clude that the ground-state of H{e = 0, A 7^ 0), which is 
always in delocalization, owns the entanglement increas- 
ing due to a < ac with ac — )■ 00, and the ground-state 
of H{e 7^ 0, A = 0) is always localized with the entangle- 
ment decreasing and in most cases with disentanglement 
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FIG. 4: (color online) The scale Oc, where (a) as a function of 
s under different characteristic parameters shows the smooth 
variance for all types of the bath; (b) as a function of s com- 
pares conventional spectral function with the modified one for 
K — 10^^^; (c) in variation with k for different types of the 
bath shows the singularity at e = 0; (d) as a function of s 
compare our displaced coherent-state approach (DCSA) for 



K = 10" 



with NRG and QMC for k = 0. 



due to the parity breaking regarding 11//. In this sense, 
any characteristic parameter calculated under e = and 
e 7^0 should be very different. So the fitting in Fig. 4(d) 
for our approach using a negligibly small k with respect 
to the NRG and QMC at k = gives the quantitative ev- 
idence that both the NRG and the QMC suffer from the 
infrared divergence for s <1 with the uncertainty equiva- 
lent to the effect of k = 10^^° in the calculation without 
the infrared divergence. 

The feature of the scaling can also be reflected in entan- 
glement. We denote the entanglement by von Neumann 
entropy E = —p^log2P+ — p-log2P~ [2^ with 



1 




which reflects the bipartite quantum correlation between 
the spin and the bath. We plot in Fig. 5 the entan- 
glement with respect to a" for different k, where af3' in 
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FIG. 5: (color online) Entanglement in variance with a" for 
different k. The curves remain unchanged for any value of s. 

In comparison with [13, [l^ for the relationship of 
the von Neumann entanglement entropy with the QPT in 
the SBM, no cusp-like behavior happens in our work for 
the entanglement changing with respect to a" and our 
results could be applied to all types of the bath. This 
is understandable because what we demonstrate is the 
scaling behavior around the critical points for the parity 
breaking, instead of the QPT for localization. Never- 
theless, similar to the results in [iJ: [3], we also find 
that the maximal entanglement appears when approach- 
ing the point a — ac from the delocalization side and 
then a rapid disentanglement at the localization side. 

In summary, we have indicated the parity breaking in 
the SBM and investigated analytically the scaling behav- 
ior of the magnetization and the entanglement as well 
as their relationship in the neighborhood of the critical 
points for the parity breaking after suppressing the intrin- 
sic infrared divergence in the spectral function. We argue 
that the conventionally employed spectral function is not 
fully reasonable and the previous conclusions drawn for 
the QPT happening in the Ohmic and sub-Ohmic SBM 
need more serious reexamination. Our analytical treat- 
ment for the scaling behavior is suitable for all types of 
the bath and should be of general interest, which is help- 
ful for clarifying different numerical results in previous 
publications and for understanding the phenomena due 
to parity breaking and the physics hidden by the infrared 
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Supplementary Material 



Eqs. (10) and (11) are in principle solvable, but time- 
and resource-consuming using currently available com- 
puting technology. For our purpose, under the condition 
A ^1, the terms of -D{m.n} with {m} ^ {n} play negligi- 
ble roles in the equations compared to other terms with 
{m} = {71} (The validity of the negligence of those terms 
is tested numerically below in Figs. 6 and 7). So Eqs. 
(10) and (11) can be reduced to 



I. ANALYTICAL SOLUTION TO SPIN-BOSON 
MODEL 



The SBM hamiltonian is given by [2[ 
a = |o-2 - -^CT:! Wfca^flfe +^ Afc(4. -f-afc)CT2, (9) 



where Oz and are usual Pauli operators, e and A are, 
respectively, the local field (also called c-number bias 0) 
and tunneling regarding the two levels of the spin, and 
ah are creation and annihilation operators of the bath 
modes with frequencies Wfc, and is the coupling be- 
tween the spin and the bath modes, which is governed 
by the spectral function J{ui) = X^fc ~ "^fc) f'^'" 
< w < with the cutoff energy 

We suppose the eigenfunction of Eq. (9) to be 



I*) 



Ew(-i)^^"'"+'rfwl{"}>j 



where C{„} and c?{„} arc coefficients to be determined 
later and \n\ = ni,--- ,nN are for different Bosonic 
modes. |{"-})a(_b) is the product of displaced coher- 
ent states of different modes [3, i.e., \{n})A{B) = 
nf=i where 



rik A, 



nk B, 



{al+qkre-'^>'''l\Q), 



o-qi/2 



with the displacement variables qk = Xk/^k and k 
1, 2, • • • , N. Using Schrodinger equation, we have 



mk-ql) + ^ 



A 



[^^k{mk-ql)' 



i{m}- 



\ ^ 

"9" «{n}£'{m,n} = -Bc{„j} , 



{»} 



(10) 



{m,n} 



Ed 



{»} 



where is given by [13, [T| 



N min[mfc,nfc] 
k=l j=0 



{m}i 



(11) 



{mk - j)l{nk - j)ljl 



k 



(12) 



(13) 

from which we may straightforwardly obtain the analyti- 
cal expressions of the cigenenergies ^{^j. = X^fe ^k(jnk — 

ql) ± ^Je^l^~^D^^~^/2, and the coefficients c^^j 

f'U/^J^ + it^Ur and = l/yiTo^ with 

/AD{ 

m,m} ■ 



II. 



VALIDITY OF THE TRUNCATION FOR 
THE BOSONIC MODES 



In the latter half of the manuscript, we calculate the 
scaling behavior of the magnetization and the entangle- 
ment by only considering the ground-state of the bosonic 
field, i.e., '^rik = 0. To check if this truncation works 
well in the case of small A, we have calculated the magne- 
tization with the bath modes of the discretizcd spectrum, 
based on the NRG logarithmic discretization |5|, where 
we used our modified spectral function and compared 
different truncations of the bosonic modes. Figs. 1 and 
2 present that the magnetization remains unchanged un- 
der different truncation of the bosonic modes for different 
bath types, which indicate our calculation using only the 
ground-state of the bosonic field to be in saturation for 
the problem. So we may consider the results based on 
our analytical treatment with ^ = to be reliable in 
the case of very small A. Moreover, in the calculations 
for Figs. 1 and 2, we employed Eqs. (12) and (13) in the 
case of X) '^fc =0, and Eqs. (10) and (11) for the case of 
J2^k =1, 2 and 3. So the good fitting of the curves in 
the figures is also the justification of the approximation 
made in Eqs. (12) and (13). 



HI. CALCULATION FOR BATH MODES OF 
THE DISCRETIZED SPECTRUM 

For the bath modes of the discretizcd spectrum, we 
modify the parameters in Eq. (9) by the NRG logarith- 
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FIG. 6: The unchanged magnetization with respect to a un- 
der sub-Ohmic dissipation s = 0.2 for difTerent truncation of 
the bosonic modes, where A = 2, luc = 1 and we have consid- 
ered the total bosons to be 0, 1, 2 and 3, respectively. 



FIG. 8: (color online) qt as the function of k under the treat- 
ment of the NRG logarithmic discretization using J(aj) (the 
curve with red dots) and J'{ui) (the curve with blue circles), 
where s = 0.1, a = 0.05, P = 10®, = 1 and A = 2. The 
inset is for s = 1 with the same values of a, P and A. 
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FIG. 7: The unchanged magnetization with respect to a un- 
der super-Ohmic dissipation s = 1.2 for different truncation 
of the bosonic modes, where A = 2, ljc = 1 and we have 
considered the total bosons to be 0, 1, 2 and 3, respectively. 



mic discretization as 



Wfe = 6 = 7fc / ujJ{io)dLO 

and Afc = 7fc/2A/7r with 7^ = Xf-cfc+^t^^ J{Lo)du. In such 

a case, it is easy to find X]/c 9fc = 7fc/(47rCfc)' which, 
as demonstrated in Fig. 8 for s < 1, is divergent with 
J{uj) but convergent using J'{uj). So our modified spec- 
tral function could effectively suppress the infrared di- 
vergence in the study of the SBM. 
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